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ABSTRACT 

We report on candidate active galactic nuclei (AGN) discovered during the monitoring of ^500 
bright (r < 18 mag) galaxies over several years with the Kepler Mission. Most of the targets were 
sampled every 30 minutes nearly continuously for a year or more. Variations of 0.001 mag and often 
less could be detected reliably. About 4.0% (19) of our random sample continuously fluctuated with 
amplitudes increasing with longer timescales, but the majority are close to the limits of detectability 
with Kepler. We discuss our techniques to mitigate the long term instrumental trends in Kepler light 
curves and our resulting structure function curves. The amplitudes of variability over four month 
periods, as seen in the structure functions and PSDs, can dramatically change for many of these 
AGN candidates. Four of the candidates have features in their Structure Functions that may indicate 
quasi-periodic behavior, although other possibilities are discussed. 

Subject headings: galaxies: active, galaxies: photometry 


1. INTRODUCTION 

A general characteristic of emission from the centers 
of galaxies that contain active galactic nuclei (AGN) is 
large amplitu de variability at all wavelengths from hour s 
to years (e.g., lUlrich et al.lll997t IWebb fc Malkanll200Q[) . 
No other phenomenon in galaxies is sufficiently powerful 
to significantly alter the brightness of an entire galaxy 


on such long timescales. Hence, surveys of , 
signed to discover AGN based on variability ( 

galaxies de- 

iHook et al.l 

1994 iBershac 

V et al. 

19981: iKlesman fc Sarajedim 2007: 

Trevese et al. 

200Hm 

acLeod et al. 20101: llVillforth et al. 

2010 c Saraiedini et al. 

2011 : Cicio et al. 20141) have been 


fairly successful. Unfortunately, ground based time series 
usually suffer from interruptions due to weather, day¬ 
light, and telescope availabilities. Ground based studies 
also suffer from atmospheric seeing which limits photo¬ 
metric accuracy, typically to several x 0.01 mag, and 
can also resul t in false detections. The Kepler satellite 
(iBorucki et al.fcOl O), which can simultaneously monitor 
a large number of galaxies at a steady 30 minute ca¬ 
dence for months without disturbance from the atmo¬ 
sphere overcomes these limitations. Detailed time series 
analysis benefits from these long continuous series. 

Although the pixels in Kepler’s CCDs subtend 4" on 
the sky, and hence the resolution is a factor 2-4 worse 
than typical ground based measurements, being in space 
compensates for this: Kepler achieves photometric accu¬ 
racy of better than 0.001 mag per 30 minute measure¬ 
ment at 17th mag and an order of magnitude better on 
timescales of ~ day. For an L* galaxy (L ~ 2.6 x 10 10 L©, 
iSchechterl ( 1976 )) with a typical central black hole of 
a few times 10 6 M 0 , the Eddington luminosity (3.2 x 
10is about as bright as the galaxy. AGN 
optical light output typically vary by large factors, and 
the ampli tude of variability anti-correlates with AGN ln - 
minosity (|Barr fc Mushotzkvlll986l: ICrisiani et al.lll996l) . 
Thus, Kepler data can reveal AGN with luminosities be¬ 
low 10 -3 of their Eddington Limit, well below values 
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reported previously. However, this technique is fairly 
insensitive to Seyfert 2 type nuclei because the regions 
immediately around their BHs are not accessible in the 
optical. 

We reduced data on ~ 500 anonymous galaxies that 
Kepler monitored during the period from Quarters 6 to 
17, a duration of about 32 months. Only a few galaxies 
were monitored throughout the entire period, but most 
have at least one year of data. In Quarter 17, Kepler 
suffered a second failure of its reaction wheels that ended 
the initial Kepler program. However, the satellite has 
been resurrected as the new K2 Mission which relies on 
the two functioning reaction wheels, but is restricted to 
pointings in the ecliptic plane. 

The optical variability of AGN is probably related to 
accretion disk instabilities, variations in accretion rates 
or changes in accretion disk structure. The emitting re¬ 
gions are generally too small to be spatially resolved, 
therefore indirect means have been found to probe ac¬ 
cretion disk dynamics and structure. Analysis of vari¬ 
ability leads to characteristic timescales that may be re¬ 
lated to light travel times across emitting objects or to 
dynamical timescales such as orbital times. Typically, 
these timescales are set by the black hole mass and the 
accretion processes (jFrank. King fc Rainell2QQ§ ). Rough 
estimates for some of these timescales are: light-crossing 

ti = 2 . 6 M 7 R 100 hr, dynamical td yn = IOM 7 -R 100 d and 
thermal cooling t t h = hr, where Mj is 

the black hole mass in units of 10 7 M©, R\oo is the ra¬ 
dius of the emission region in units of 100 times the 
Schwarzschild radius (2GM/ c 2 ), and v is t he Shakura- 
Sunyaev viscosity parameter (jShakura fc Sunvaevlll973[ ). 
For typical AGN models, these timescales range from 
hours to years. It is expected that, with sufficiently com¬ 
plete and accurate measurements, fits to the Power Spec¬ 
tral Density (PSD) function or Structure Function (SF) 
of an AGN’s light curve might display features at these 
key physical timescales. 

Several teams hav e published results from Kepler 
data of known AGN (jCarini fc Rvlell2012t Chen fc Wand 
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Fig. 1.— Distribution of the galaxy sample on the sky in equato¬ 
rial coordinates. The symbols tell in which years the galaxies were 
observed. Cycle 2, squares; Cycle 3, crosses; and Cycle 4, pluses. 
The points are overlying a Palomar Sky Survey image of the Kepler 
Field. 

120151 : iRevalski et al.l I2014D . iMushotzkv et al.l (|2011l) re¬ 
ported on four AGN, finding a powe r-law slo pe for the 
PSD with steep (-2.6 to -3.3) slopes. IWehrle et al.l (|2Q13f ) 
have found PSD function fits by power laws with no clear 
breaks for four other AG N. Th e PS D o f acti ve galaxy 
Zw 229-15 was studied bv lEdelson et al.l (j2014f h and was 
found to have a strong bend at the ~ 5 day timescale, 
which is similar to what has been reported in the x-ray 
PSD of AGN of similar mass. 

Here, we try to make use of Kepler’s exquisite photo¬ 
metric capabilities to explore nuclei with activity levels 
lower than would normally be detectable and to focus on 
determining the true occurrence rates of AGN in galax¬ 
ies. In §2, we describe the target selection procedure 
and give some characteristics of the sample of galaxies 
observed. In §3, our photometry procedure is discussed. 
The long term instrumental trends are removed from our 
aperture photometry by making use of the project sup¬ 
plied co-trending basis vectors. In §4, a method to dis¬ 
criminate for activity is described and applied to find 
candidate AGN and structure functions are made from 
those light curves and analyzed. 

2. OBSERVATIONS 
2.1. Source Selection 

The Kepler Field, centered at galactic l = 264?5, b = 
+ 13? 5, lies at fairly low galactic longitude and therefore 
had not been well covered by AGN surveys. Pre-launch 
of Kepler, only a few AGN were known that resided in 
the field. Soon after launch, there were efforts to search 
surveys at various wavelengths to identify AGN based 
on colors (lEdelson fe Malkanll2012f ) or x-ray flux. But, 
for our purpose, to discover the true occurrence rate of 
AGN, we simply selected a set of undistinguished galax¬ 
ies of small angular extent, mostly from the 2MASS sur¬ 
vey, to be monitored with Kepler. Because the data 
rate from the Kepler camera was so large and the down¬ 
load bandwidth was limited, specific targets in the Ke¬ 
pler Field were competitively awarded each year, and 
only a small patch of CCD pixels around each target 
image was downloaded. Therefore the galaxies had to 
be small in angular extent to minimize the number of 
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Fig. 2. — The number of galaxies monitored per Kepler Quarter 
and assigned to this team’s programs. 



Number of Kepler Quarters of Monitoring 


Fig. 3. — Histogram showing how the target galaxies are dis¬ 
tributed by the number of Kepler Quarters they were observed. 


pixels for downloads. Our target galaxies (in programs 
G020058, G030032, and G040057) were e xtracte d from 
the 2MASS extended source catalog (2MX; iJarrett et al.l 
( 2000 )). The target apparent Kepler magnitudes were 
restricted to < 17.3 mag, and we selected small galaxies 
as defined by the 20th mag K-band isophotal radii, in¬ 
cluding 3cr errors, 95% were < 2 pixels and the average 
was 3.7 ±1.3 arcsec, or 0.93 ± 0.33 Kepler pixels. Ap¬ 
proximately every 3 months the spacecraft was rotated 
90° to maintain sunlight on the solar panel and to down¬ 
load data. This defines a Kepler Quarter. Targets were 
deselected if they fell between CCDs after any rotation. 
According to the Explanatory Supplement (Cutri et al 
2000), at these low galactic latitudes, about 80% of the 
extended 2MX sources are galaxies and most of the rest 
of the sources are multiple stars. Our selection of 2MX 
galaxies contained only sources that were visually con¬ 
firmed by the 2MX authors. This procedure resulted in 
1,858 galaxies as bright as 10th mag, and the median was 
near 15 mag. We then selected 400 galaxies from this set 
that were best separated from other sources in the cata¬ 
log and least contaminated from the overlapping scaled 
point spread functions of neighbors. Their distribution 
on the sky is show in Fig. [lj 

Our list of galaxies was proposed in Cycle 2, but only 
130 were accepted for downloading and some for less than 
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4 Quarters. In Cycles 3 and 4 most of the galaxies on 
the list were selected for downloading. Figure [2] shows the 
number of galaxies awarded to this project and actually 
monitored in each Quarter. For two years about 400 
galaxies were monitored, and in total, 519 galaxies were 
assigned to be monitored. Quarter 17 was cut short due 
to the failure of a reaction wheel and the photometry 
in that quarter appears to have suffered instabilities in 
most cases. Figure [3] presents the number distribution 
of the number of Quarters observed. The distribution of 
the sample on the sky is shown in Figure [U 

3. DATA REDUCTION 

The Kepler science data pipeline begins with a cal¬ 
ibration (perf ormed by the CAL module) of the CCD 
images (jQuintana et all [2010) by correcting for; bias 
level, dark current (insignificant), smear due to shutter¬ 
less readouts, nonlinear gain (median of 112 e - /ADU), 
certain undershooting pixels, and individual pixel sen¬ 
sitivities (flat field). The pipeline also produces cal¬ 
ibrated light curves in its Photometric Analysis (PA) 
package, based on the total counts above background 
in an ’’optimal aperture” around the target defined by 
the project. The optimal aperture is a set of contiguous 
pixels within the downloaded stamp, fixed throughout a 
Quarter, with count rates from the target high compared 
to the background but avoids possible light contamina¬ 
tion from nearby sources. Adjustments are made for cos¬ 
mic rays and the observation time array is converted to 
barycentric times. Arrays of centroid positions and back¬ 
ground light as functions of time are also created. The 
next sta ge of red uc tion , the Presearch Data Conditioning 
(PDC) (|Stumpe et al.l 20121 ). applies corrections based 
on both known instrumental and spacecraft anomalies 
and correlated systematic photometric excursions found 
in the data. Excess flux in the source apertures due to 
unavoidable contamination from nearby stars is removed 
and outlier data points are flagged. 

Variations in PA counts are found to correlate with fo¬ 
cus and PSF changes which arise from flexing of the tele¬ 
scope structure. Thermal transients, lasting several days, 
occur following downlinks, after safe modes, or from the 
changing orientation of the spacecraft with respect to the 
Sun. Drift in pointing at the sub-pixel level also mod¬ 
ifies the counts. Differential velocity aberration across 
the large field unavoidably causes the position of target 
images to slowly change, especially at the edges of the 
field. In the PDC process, a set of 14 normalized vectors, 
called cotrending basis vectors (CBV), are generated by 
singular value decomposition that represent correlated 
instrumental artifacts in the LCs in each CCD and in 
each Quarter. For the PDC light curves, a superposi¬ 
tion of these vectors that minimizes the rms deviations 
for a target over the Quarter is subtracted from the light 
curves. 

It was recognized that the procedure of minimizing 
deviations could inadvertently remove astrophysical fea¬ 
tures in the time series and therefore, whenever the vari¬ 
ations exceed some strength and are clearly intrinsic to 
the target, this step is either skipped or fewer CBVs are 
used. For known AGN (in other GO data) with large am¬ 
plitude (> 5%) variations, the PDC light curves supplied 
by the Kepler Project seem reasonable, and it appears 
that either no CBVs, or just the first one or two, are 


used. However, in some of our target galaxies with lower 
amplitude variations, the PDC procedure may have re¬ 
moved some of the long timescale physical fluctuations. 
Usually, one can still discern that these are active galax¬ 
ies because the rms in the light curves remains elevated 
compared to the very quiet galaxies, and power in the 
PDC is shifted from longer to shorter timescales. Other 
examples of PDC over correction are the super novae dis¬ 
covered in our programs (Oiling et al 2014) in which the 
PDC light curves are so distorted that the supernova 
events were rendered nearly unrecognizable. 

Another problem is that many quiet galaxies have light 
curves that appear active in certain Quarters. In a few 
cases, this occurs because the PDC procedure just failed. 
In other cases, the CCD had undergone a period of excess 
noise arising from crosstalk between the CCDs and/or 
the 4 fine guidance sensors or from a high frequency 
amplifier oscillatio n as has been documented for Kepler 
CCDs (jKolodzieiczak et al.lll2010 ). In most of these cases 
it is necessary to simply disregard entire Quarters for tar¬ 
gets on these CCDs to ensure reliable measurements. 

3.1. Generating Long Term Light Curves 

We reprocess our galaxy data into light curves starting 
from the Kepler pipeline CAL output of postage stamp 
cutouts with the background subtracted. By using more 
pixels than the PDC procedure, the long timescale trend¬ 
ing effects from pointing drift and focus changes are re¬ 
duced at the expense of some additional readout and sky 
noise. For our purposes, this tradeoff is highly beneficial 
because the former greatly exceeds the latter. Therefore, 
we sum counts in either 3 x 3 or 5 x 5 pixel apertures. 
With 5 x 5 pixels, the trends over a Quarter, which typi¬ 
cally have an amplitude of 3 - 4% in the PA light curves, 
drop to under 1%. We do not use 7x7 pixel apertures 
because they often suffer significantly from contaminat¬ 
ing variations of neighbors and from larger errors arising 
from a given background error. 

We looked at two methods for correcting the long term 
trends in the counts due to instrumental effects. One 
method is to determine PSF widths and centroids in two 
dimensions and then find first order correction factors 
on the four parameters. We could have used the Ke¬ 
pler pipeline determined widths and centroids, but it is 
easy to calculate simple stand-ins for these measures. For 
width changes we use the ratio of counts in the 5x5 pixel 
aperture to that in the 3x3 pixel aperture, and for cen¬ 
troid shifts in each direction, we use the difference in the 
counts in the pixels on either side of the center divided 
by the three contiguous central pixels. Since the trends 
in counts have amplitudes of only a few percent, it is rea¬ 
sonable to assume that terms linear in these parameters 
would dominate. By assuming the coefficients are con¬ 
stant over each quarter, one can simultaneously solve for 
coefficients that minimize the variations. For the quiet 
galaxies in our sample, this procedure typically reduces 
the trends to less than 0.1% variability over 3 months. 
This confirms that much of the trending in Kepler pho¬ 
tometry is due to focus changes and centroid motion as 
described above. 

However, we found a better procedure with which we 
base our conclusions. We apply the project’s CBVs to 
our aperture photometry, using an altered goodness cri¬ 
terion for the coefficients to create a LC with normalized 
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corrected counts (NCC): 


NCC = 


Counts 

(Counts) 


- ZiCi x CBVi 


A criterion based on minimum rms deviation from the 
mean can falsely accept a solution in which the corrected 
light curve oscillates around the mean more, yet fewer os¬ 
cillations is preferred. To understand this better, imagine 
a signal that rises near the beginning and then becomes 
constant and also imagine a single CBV that is just an 
upwards convex function, ie. does not match the shape 
at all. Subtracting a multiple of the CBV could reduce 
the rms by a substantial amount by adding two more 
crossings of the average, but since it does not match the 
pattern of variability, it should be given a very small co¬ 
efficient. What is required is a criterion that increases 
the coefficient if doing so improves the flatness of the 
corrected signal everywhere locally in time and not by 
having more zero crossings. We find good stable results 
with a criterion in which we minimized on the sum of 
the squares of the changes separated by any length from 
10 - 25 days. If a shorter delay is used then the low 
frequency trends are not sufficiently diminished, and if 
a longer delay is used then the low frequency trends are 
over suppressed at the expense of the high frequencies. 
We also limit the value of the coefficients that we apply to 
the normalized CBVs to < 1.5 to prevent over-correction 
on the normalized LCs. In the PDC procedure, these co¬ 
efficients have been found to typically be under 1.0. We 
start by using the first two CBVs and include an addi¬ 
tional CBV if the decrease in the criterion was significant, 
up until a maximum of six CBVs. The Quarterly CBVs 
provided by the project are ordered by the amount of 
correlation that they have with all of the LCs of targets 
on a given CCD of the camera. 

After the correction are made to a Quarter, we sew 
pairs of Quarters together by fitting a few days at each 
Quarter’s border with third degree polynomials. A time 
is selected several days into the later Quarter and the 
ratio of the two polynomials at that time is the constant 
multiplier applied to the later Quarter. 

After about one year, Kepler returns to the same Sun- 
Kepler-FOV configuration and the counts and the trend¬ 
ing pattern repeats fairly closely that of the previous 
year. This can be useful for solving the long term trend¬ 
ing puzzle. We first use the 3x3 pixel apertures, and 
where we have observations over more than one year, we 
check that the fluxes now repeat well even without sewing 
every 5th Quarter to its previous one. If the yearly repeat 
is off by more than the error of a single exposure, then 
we redo the photometry with a 5 x 5 pixel aperture, and 
this usually works well. If there is still a problem, often 
it is an indication that we have overused the CBVs and, 
usually, reducing the number of CBVs utilized achieves 
good matches on the one year timescale. 


4. NEWLY DISCOVERED ACTIVE GALAXIES 

After correcting for long term trends and eliminating 
Quarters where CCD bias artifacts were dominant, only 
4 active galaxies were easily discernible by visual inspec¬ 
tion of the all the light curves (seen in Fig. [5] and [6j). To 
detect lower amplitude active galaxies, we made use of 
the fact that for known AGN, typically, the PSDs, fit on 



Fig. 4.— Ratio of Structure Function at two delay times, 15 
days and 5 hours, as a function of Structure Function at 5 hours 
for 457 galaxies measured by Kepler Mission. The candidate AGN 
are numbered 1 to 19 (identified in Table 1) and have ratios > 3x 
the standard deviation of the distribution (dashed lines), cr = 0.11. 
The most active galaxy lies off the plot at a ratio of 29.9. 

the timescale of weeks or longer, have power laws going 
as f~ a with a between 1.5 and 2.5. Thus, the ratio of 
the power at many days to that at just a few hours could 
be a discriminator for detecting AGN. 

To obtain a statistically valid estimate of the active 
galaxy occurrence rate, we devised an automated pro¬ 
cedure that would resist being triggered by temporary 
periods of high instrumental noise occurring in many of 
the CCDs. The sample was restricted to the subset of 
targets with three or more Quarters of observations; 473 
galaxies remained in the set. For each of these, struc¬ 
ture functions were made on the quietest Quarter of the 
PDC LC; this ensures that none are being selected as 
active candidates because of excess instrumental noise. 
Although this precaution may result in excluding active 
galaxies that happened to be very quiet during one 3- 
month period, we believe that this is rare, while the dan¬ 
ger of being fooled by instrumental noise was high. We 
defined a dimensionless, unweighted structure function, 

([F(* + T)-fV)] 2 
' ( Fit )> 2 

where F(t) is the brightness at time t, r is a delay time, 
and angle brackets mean average over all points in the 
time series. We first smooth the data by averaging sets 
of four 30 minute exposures, in this form, the structure 
function directly provides the mean square of the frac¬ 
tional change at different delay times in 2 hour flux mea¬ 
surements. The angle brackets indicate an average is 
taken over the differences in all pairs in the LC at a 
given delay time. 

For galaxies with no activity, the SF should be a con¬ 
stant at all delay times. As an automated procedure to 
discriminate active from quiet galaxies, we form the ratio 
of the SF with delay times at 5 hours to that at 15 days, 
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TABLE 1 

Active Galaxies Discovered in sample of 474 galaxies 


Rank 

Kepler ID 

RA 

Dec 

g 

r 

J 

Skygroup 

SF(10d)/SF(4h) 

1 

10645722 

18 47 22.340 

+47 56 16.13 

15.958 

15.635 

15.505 

11 

29.95 

2 

5686822 

18 59 12.137 

+40 59 26.73 

16.801 

15.872 

15.520 

50 

3.66 

3 

11716536 

19 35 17.046 

+49 51 51.67 

17.063 

16.218 

16.302 

14 

3.35 

4 

7986325 

19 57 55.472 

+43 46 49.27 

16.850 


15.456 

53 

2.54 

5 

6751969 

18 47 30.765 

+42 17 07.23 

16.002 

15.515 

15.594 

31 

1.99 

6 

7691427 

19 38 47.414 

+43 21 57.31 


17.114 

15.080 

58 

1.85 

7 

12553112 

19 12 52.522 

+51 35 06.37 

17.844 

16.666 

15.831 

3 

1.66 

8 

11768473 

19 38 04.715 

+49 55 47.89 

16.222 

15.722 

15.936 

14 

1.56 

9 

2142191 

19 06 03.085 

+37 35 09.70 

15.696 

14.759 

15.021 

71 

1.52 

10 

12556836 

19 22 00.693 

+51 35 56.09 

14.785 

14.321 

14.939 

4 

1.51 

11 

11808151 

19 14 16.560 

+50 00 24.84 

18.229 

17.251 

16.356 

2 

1.50 

12 

5683305 

18 52 03.013 

+40 59 33.73 

17.005 

16.274 

15.714 

52 

1.47 

13 

8024526 

19 21 25.589 

+43 53 16.90 

14.971 

14.393 

14.282 

41 

1.45 

14 

5511084 

18 51 55.133 

+40 47 23.88 

14.984 

13.925 

14.442 

52 

1.44 

15 

6714622 

19 50 09.928 

+42 06 53.09 

17.544 

16.228 

15.064 

74 

1.42 

16 

4148802 

19 16 58.775 

+39 14 48.41 

18.802 

17.757 

15.789 

68 

1.40 

17 

10402746 

19 24 07.041 

+47 33 41.50 

16.991 

16.245 

15.855 

18 

1.40 

18 

9509125 

18 50 29.335 

+46 08 01.14 

18.376 

17.357 

15.487 

10 

1.40 

19 

8884097 

19 25 08.969 

+45 10 29.57 


17.699 

15.601 

43 

1.38 


averaging every four time samples to reduce noise. This 
short delay time is almost as short as we can go with 
the averaging used. At these particular two points in the 
delay times, the number of standard deviations from the 
mean in this ratio are maximized for these candidates 
and subsequently the maximum number of candidates is 
found. Figured] shows the distribution of this ratio as a 
function of the SF at 5 hours. The mean of the ratio was 
1.03 with root mean square deviation of a = 0.11. Of 474 
galaxies, 19 galaxies have ratios > 3cr, where a is the root 
mean square deviations of the ratios with 5a deviants re¬ 
moved, and these are labeled in the diagram with their 
rank order. In a normal distribution, only 0.15% of the 
sample would be greater than 3a from the mean, which 
for this size sample corresponds to 0.7 galaxies. 

Presently we do not know much about these AGN can¬ 
didates. Searches on NED and HE AS ARC show that all 
but two of the 19 candidate AGN are unstudied other 
than by 2MASS measurements. None have detected 
x-ray emission. KIC 12556836, also known as NVSS 
J192200+513557, was observed to emit 4.2 ± 0.5 mJy at 

1.4 G Hz in the NRAO VLA Sky Survey (jCondon et al.l 
Il998[h but no spectra have been taken. KIC 8024526 
was determined to have a redshift of z = 0.0551 (m - M 
= 36.80, Ho = 72 km/s/Mpc) ([Faber fe Dressier! 1 19771 ) 
and thereby confirmed to be a member of Abell Cluster 
2319. The absolute magnitude for the whole galaxy is 
M r = -22.4 mag (3.8 x 10 10 L o = 2.2 x 10 44 ergs/s ), 
and so the detection of variability at the level of 0.5 
milli-magnitudes corresponds to luminosity variations of 
1.1 x 10 41 ergs/s at 230 Mpc. 

KIC 10645722 (the most active candidate, Fig. [5]) was 
selected to be studied spectroscopically at Lick on the 
basis of its variability on Kepler and having infrared col¬ 
ors from WISE in the region pred ominantly populated 
by AGN (lEdelson fe Malkanl 1201 2). There is an r-band 
CFHT Megacam image of it, and it was imaged in Cycle 
21 on HST with the WFC/UVIS (HST proposal 13409) 

. It is a Seyfert 1 at redshift z=0.068 (K.L. Smith, per¬ 
sonal communication) or distance modulus 37.26 (Ho = 
72 km/s/Mpc). This galaxy has a total luminosity of 

1.5 x 10 44 ergs/s and HST and MegaCam images show 


that the majority of its light comes from the central 
source. 

Variations are seen, over three years, having amplitude 
of ± 4% for this source, while other candidates vary by 
< ± 0.5%. Assuming the BH is emitting below the Ed¬ 
dington Limit (LEdd = 1.3 x 10 38 Mbh/M q erg/s), the 
minimum mass of the BH is ~ 1.1 x 1O 6 M 0 . 

If we identify elliptical galaxies as those having photo¬ 
metric color r — J > 1.0 mag, then four of the 19 are el¬ 
lipticals: KIC 7691427, 6714622, 4148802, 9509125, and 
8884097. Interestingly, these are all toward the lower end 
of the ranking in the detection ratio with 3 of the 5 in 
the bottom 4. Given the favorable weighting of elliptical 
galaxies in the selection criterion for our targets, this in¬ 
dicates a strong preference for detecting spiral galaxies 
as active. 

Unfortunately, we do not have redshifts yet for most 
of these galaxies, so we cannot convert their brightnesses 
to luminosities. 

4.1. Structure Functions and PSDs 

To characterize the variations as a function of 
timescale, we use the unweighted Structure Function de¬ 
fined above. SFs are stable if the time series contains 
many gaps or long gaps, as long as the data are not too 
sparse. Small gaps exist in this data from several causes, 
chiefly: cosmic ray events, several day events to turn the 
spacecraft and download the data to earth on a monthly 
basis, Quarterly rotations of the spacecraft, Moire pat¬ 
tern noise and other electric interference events on the 
CCD. Sometimes a given target was simply not sched¬ 
uled for a downlink for a specific Quarter. 

The top two panels of Fig. [5] - [14] show LCs of rep¬ 
resentative Quarters for each candidate. The LCs, af¬ 
ter our processing (black lines) differ by varying degrees 
from the Kepler project’s processed LCs (red lines). In 
some cases, the PDC failed completely or nearly so in 
the face of large variability. The SFs for all quarters 
with acceptable quality data are plotted (black for our 
reduction, red for the project’s PDC, and blue and vi¬ 
olet for the Quarter of the LC presented) in the lowest 
panel of these figures. However, if the SF of the PDC 
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LC is unreasonably noisy compared to our LC, it is not 
plotted. The SFs are normalized to have the same am¬ 
plitude at the highest frequency. These galaxies all show 
excess noise at delay times of about a day that further in¬ 
crease for longer delays, as is typical of previously known 
AGN. However, from one Quarter to another the overall 
amplitude of the SF can change significantly. Thus the 
activity in these cases is not consistent with the simple 
auto-regressive or damped random walk models. This 
was found by iKasliwal et al.l ( 2015 ) to be the case for 
the more dominant AGN in Kepler studies. 

In most of the SFs, the power at the shortest delay 
times measured (an hour) is dominated by the Poisson 
noise of the target and sky and instrumental noise such as 
readout noise, and this explains the tendency to flatten 
out at delays of a day or less. The brightest sources, 
though, clearly show physical activity registering all the 
way down to one hour. In most of the cases, a knee 
occurs in the SF between 2 and 20 days where the LC 
nearly flattens to zero slope, indicating little additional 
variability at long timescales. 

The second most active, KIC 5686822, has SF power 
law index of ~ 1.0 around 1-6 day delay times, yet 
the index falls to nearly zero by 10 day timescales. LC 
excursions over ^ 6 days are typically at the 5 mmag 
level. 

Structure Functions react to periodic behavior by drop¬ 
ping to lower values at integer multiples of the period. 
Several of the SFs consistently show clear and deep min¬ 
ima (KIC 10645722, KIC 5686822, KIC 7986325, and 
KIC 4148802), but the amplitudes and delay times of 
the minima change over the observing. An exception is 
KIC 7986325 which appears to maintain a single period 
but has big amplitude swings. It is conceivable that this 
one is contaminated by a variable star. We are working 
on simulations to determine if this quasiperiodic behav¬ 
ior is likely to be statistical artifacts arising from power 
at longer frequencies (red leak), or from residuals of the 
trends in the Kepler data, or real quasiperiodic behavior 
in the accretion disks of BHs on the time scales of 20 - 40 
days. In the case they are Kepler artifacts, these galax¬ 
ies remain AGN candidates since astrophysical activity 
is the best explanation for why these LCs could not be 
fully detrended by the cotrending analyses even though 
the amplitudes of variation are relatively small. 

The conversion from SF to PSD is complex 
unless they are both simple power laws, there¬ 
fore we have also created PSDs (online data, 
http://www.astro.umd.edu/^eshaya/AGN/) from our 
LCs. These also flatten at frequencies above v = 1/day 
or so, and rise to lower frequencies, but most begin to 
flatten at 10d -1 . A table (also online) gives the param¬ 
eters of fits {yb r k , /3 ) to the broken power law: 

PSD(v) = PSD^Vbrk) ^ for o.ld -1 < V < V^rk 

\ Vbrk J 

( V \-P 

= PSD(u brk ) - for v > v brk 

\ Vbrk J 

where is the frequency of the break from the white 
noise regime to the target’s variability signal. At fre¬ 
quencies below O.ld -1 , these PSDs are often not fit well 
by a power law. The fit parameter values vary widely in 


many of the candidates again indicating that the ampli¬ 
tude of the variations changes over time. Of course, this 
would be expected if there is significant power on longer 
timescales. 

5. CONCLUSIONS 

AGN selection criteria have traditionally been based 
on imaging and spectr oscopic techn iques in the various 
wavelength bands (see iMushotzkvi ( 200 4) for a discus¬ 
sion). However when the light from the AGN is signif¬ 
icantly less than that from the host galaxy the unique 
AGN spectral or imaging signatures can be difficult to 
recognize and result in severe incompleteness at even 
moderate redshift. Even in the x-ray band, where the 
contrast with stellar emission processes is the greatest, 
at low luminosities (< 10 42 ergs/sec), AGN emission can 
be confused with the emission from ultra-luminous x- 
ray (ULX) sources, and at even lower luminosities with 
emission from bright x-ray binaries. The virtual ubiq¬ 
uity of optical variability from AGN combined with the 
very high signal to noise in the Kepler data open up a 
totally new channel for detecting AGN and allow direct 
comparison with other survey techniques. 

Using the Kepler Mission instrument, we measured un¬ 
precedented detailed light curves of ^500 galaxies. We 
were able to improve on the long term trending, in many 
cases, over the light curves generated by the project 
pipeline by using larger apertures and re-estimating opti¬ 
mal coefficients for the cotrending basis vector detrend¬ 
ing step. Although non-optimal for the noise level of 
the 30 minute integrations, the larger apertures (typi¬ 
cally 5x5 pixels) greatly reduced the longer timescale 
noise in the light curves by being less sensitive to PSF 
shape changes and centroid motions. The project’s de¬ 
rived cotrending basis vectors nicely removed most of the 
remaining instrumental trends. 

Of the 474 galaxies with long baselines and imaged on 
quiet CCDs, 19 galaxies (~4%) show excess activity in¬ 
dicative of having an active galactic nucleus. Specifically, 
they had more variability at 15 days than at 5 hrs. The 
other roughly 460 galaxies were quiet, varying typically 
at levels below 5 milli-magnitudes over 15 days. Many 
of our candidates have variability at levels that are not 
discernible by just viewing the LC, as our technique of 
using the SF ratio can detect quite low levels of variabil¬ 
ity. In the SFs of about half of the active galaxies, 10, 
show a knee at delay times of between 4 to 20 days above 
which the rate of rise in the amplitude of variability with 
delay times slows considerably. But for four of these, the 
cause of the decrease toward longer time scales in the 
SF is due to periodicity. It will require further analysis 
to determine if this is due to residual Kepler sensitivity 
trends, statistical fluctuations from breaking up the time 
series into short, roughly 90 day, sections, or quasiperi¬ 
odic behavior of the AGNs. 

We plan on obtaining optical spectra of the Kepler 
AGN candidates (some are already acquired) to deter¬ 
mine whether they show spectral AGN signatures, but 
even that is not a certain measure, since a fair fraction of 
x-ray selected AGN do not show eit her optical or near- 
IR AGN associated emission lines (jSmith et al.l 2014 ). 
Spectra could reveal if the narrow- or broad-line regions 
of these candidates are different in composition or struc¬ 
ture from AGN with more pronounced variability. 
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The techniques developed in this paper have great ap¬ 
plications to observations with LSST or the Zwicky Tran¬ 
sient Factory which both should have sufficient precision, 
cadence and time coverage to reveal similar variability 
characteristics to what we have found with Kepler. We 
have a program to obtain more AGN with the K2 mis¬ 
sion. While the K2 observations consist of single 70 - 
80 day time series for each field, the structure functions 
from Kepler show that this is sufficiently long to detect 
more candidates and obtain unique information on their 
power spectrum or structure functions. 

The authors were supported by Kepler guest observing 
grants NNX12AC95G and NNX13AC27G from NASA. 
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Fig. 5.— Corrected and normalized single quarter light curves (black dots) of galaxies KIC 10645722 using a 5 X 5 pixel aperture (left) and the first two CBVs and KIC 5686822 
(right) using a 3 X 3 pixel aperture and CBVs 1-3. The red dots are the project’s PDC corrected light curves using ’optimal apertures’. These two have strong enough activity to be 
obvious AGN candidates. The former needed just two CBVs for long term detrending because the amplitude of its variation was high. Below the light curves are structure functions 
from the various quarters where the target was observed and the detector was well behaved. For the quarter shown above, we use blue lines for our analysis and violet for the PDC, 
but for other quarters we use black and red. 
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Fig. 6. — Corrected light curves (top) and Structure Functions (bottom) from our analysis (blue and black) and the project’s PDC (red). KIC 11716536 (left) and KIC 7986325 
(right). 
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Fig. 7.— Corrected light curves (top) and Structure Functions (bottom) from our analysis (blue and black) and the project’s PDC (red). KIC 6751969 (left) and KIC 7691427 
(right). 
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Fig. 8. — Corrected light curves (top) and Structure Functions (bottom) from our analysis (blue and black) and the project’s PDC (red). KIC 12553112 (left) and KIC 11768473 
(right). 
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Fig. 9. — Corrected light curves (top) and Structure Functions (bottom) from our analysis (blue and black) and project’s PDC (red). KIC 2142191 (left) and KIC 12556836 (right). 
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Fig. 10. — Corrected light curves (top) and Structure Functions (bottom) from our analysis (blue and black) and project’s PDC (red). KIC 11808151 (left) and KIC 5683305 (right). 
Our reduced light curve of KIC 11808151 has less growth in fluctuations at longer delays the the PDC curve. Therefore, this is not a likely candidate. 


CO 


AGN Discovered in the Kepler Mission 






















Shaya, Oiling & Mushotzky 


KIC 8024526, 5x5 pix, CBV=1- 6 

1— 1 1 1 



0.996 =* —I * I I ■ I—I—I— J—\— h—I —I— I " l ' I— I- 
1.004 y . . 1200 . 1220 . 1240 1260 


H —H= 


» * • • 

=j_i_i_1_i_i_i_1_i_ 

_i_1_ 

J_L 

_l_ 

J 

_l_1_ ± 

1200 1220 1240 

Time [MJD - 54,833] 

- 1 —i i 111 in- 1 —i i i 111 ii- 1 —i i 111111 

1260 

— E —1-1 1 1 1 1 1 1 


CD 

i— 

=3 

-i—■ 

O 

=3 

L_ 

0) 


2.56x10 



0.01 


0.10 


1.00 
t [Days] 


10.00 


100.0 


KIC 5511084, 3x3 pix, CBV=1-6 



1480 


1500 1520 1540 

Time [MJD - 54,833] 


o 

i— 

o 

d 

i— 

if) 



% [Days] 


Fig. 11. — Corrected light curves (top) and Structure Functions (bottom) for our analysis (blue and black) and project’s PDC (red). KIC 8024526 (left) and KIC 5511084 (right). 
KIC8024526 had a glitch in Q9, the quietest rms Quarter. However, Quarter 11 is quiet, which indicates that this galaxy is not active. 
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Fig. 12.— Corrected light curves (top) and Structure Functions (bottom) for our analysis (blue and black) and project’s PDC (red). KIC 6714622 (left) and KIC 4148802 (right)The 
variability in light curves of both KIC4148802 and KIC6714622 were too suppressed by the pipeline processing. For KIC6714622 we were able to make use of one year repeatability. 
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Fig. 13.— Corrected light curves (top) and Structure Functions (bottom) for our analysis (blue and black) and project’s PDC (red). KIC 10402746 (left) and KIC 9509125 (right) 
KIC10402746 is inactive in Quarter 16. 
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Fig. 14. — Corrected light curves (top) and Structure Functions (bottom) for our analysis (blue and black) and project’s PDC (red) for KIC 8884097. 
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